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Abstract 

We study the two-dimensional Eguchi-Kawai model as a toy model of the IIB matrix 
model, which has been recently proposed as a nonperturbative definition of the type 
IIB superstring theory. While the planar limit of the model is known to reproduce the 
two-dimensional Yang-Mills theory, we find through Monte Carlo simulation that the 
model allows a different large N limit, which can be considered as the double scaling 
limit in matrix models. 
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1 Introduction 



It is no doubt that the nonperturbative study of string theories is one of the most exciting 
topic in particle physics, since it might provide all the answers to the fundamental ques- 
tions concerning our world, such as the space-time dimension, the gauge group, the three 
generations of the matter fields and so on. Indeed the nonperturbative study of the bosonic 
string theory in less than one dimension was successfully done through the double scaling 
limit of the matrix model [|l[] some time ago. Much effort has been made towards the gener- 
alization of the approach to bosonic strings in more than one dimensions |3| as well as to 
superstrings 0. f|. 

Recently, a new type of matrix model has been proposed as a nonperturbative definition 
of the superstring theory. The model proposed by Banks, Fischler, Shenker and Susskind 
0, which is conjectured to provide a nonperturbative definition of the M-theory in 
the infinite-momentum frame, is a matrix quantum mechanics, which can be obtained by 
dimensional reduction of the 10D large N super Yang-Mills theory to one dimension, while 
the one proposed by Ishibashi, Kawai, Kitazawa and Tsuchiya |§, which is conjectured to 
provide a nonperturbative definition of type IIB superstring theory, is a matrix model, which 
can be obtained by dimensional reduction of the 10D large N super Yang-Mills theory to a 
point. It is now conventional to call the former as M(-atrix) theory and the latter as IIB 
matrix model. Although these proposals are basically based on the spirit to describe the 
strings in terms of matrices of very large size, there are many novel features that didn't show 
up in the old-fashioned matrix models. One of them is the appearance of the Yang-Mills 
theory. This is related to the effective action |9j of D-particles and D-instantons, respectively 
for the two models. It can also be viewed to be related to the supermembrane action and the 
Green-Schwartz action in the Schild gauge, respectively. Another interesting feature of these 
models is that all the target-space coordinates come out as the eigenvalues of the matrices. 

For the IIB matrix model, the light-cone string field Hamiltonian of the type IIB super- 
string theory has been reproduced [JTO] by identifying the string field for fundamental strings 
with the Wilson loop and analyzing the Schwinger-Dyson equation. As an important conse- 
quence of this study, they have clarified how one should take the double scaling limit in this 
model. The IIB matrix model can be considered as an example of large N reduced models 
(See Ref. [|11|] for a review.), which have been studied so far exclusively in the planar limit, 
in which the models are equivalent to the field theory before being reduced. Whether a large 
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TV reduced model allows any sensible double scaling limit is itself a nontrivial question. 
In this paper, we study the two-dimensional Eguchi-Kawai model [ I2"|| , as a toy model of 



the IIB matrix model. The model is nothing but the SU(iV) lattice gauge theory on a 1 x 1 
lattice with the periodic boundary condition and it is equivalent to the SU(iV) lattice gauge 
theory on the infinite lattice in the planar limit. We perform a Monte Carlo simulation of 
the model and find that the model indeed allows a different large N limit, which can be 
considered as the double scaling limit in matrix models. 

The paper is organized as follows. In section 2, we give the definition of the model and 
review some known results. We also explain the algorithm we use in our simulation. In 
section 3, we study the planar limit of the model. We show how the finite N effects appear 
by measuring Wilson loops. This gives us an important clue to a possible double scaling 
limit. In section 4, we present the data which show the existence of the double scaling limit 
in the model. Section 5 is devoted to conclusions and discussions. 

2 The model 



The Eguchi-Kawai (EK) model is defined by the following action [12 



Sek = ~N(3 J2 triUpUvUpl}, (2.1) 
This model has the U{1) D symmetry. 

U, - e^U,. (2.2) 



In Ref. [|12|| , it has been shown that if the U(l) symmetry is not spontaneously broken, the 
model is equivalent to the SU(iV) lattice gauge theory on the infinite lattice in the large N 
limit, where the coupling constant (3 in the action Q2.1|) is kept fixed. This limit is referred 
to as the planar limit, since in this limit Feynman diagrams with planar topology dominate. 
By equivalence, we mean that the expectation value of the Wilson loop 

W(C) = ^tr {U x , a U x+&>/3 U x+&+ ^ ■ ■ ■ U x _ kx } , (2.3) 

which is the product of link variables along a closed loop C, calculated within the lattice 
gauge theory with the action 

Slgt = -NPY1 E tr{^[/ I+ ^^ ifl ^}, (2.4) 
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is equal to the expectation value of the corresponding observable 

w(C) = Itr {V a U p U 1 ■■■U x }, (2.5) 

calculated within the EK model with the action (|2.1|) in the planar limit. 

It is found in Ref. |y| that when the space-time dimension D is larger than two, the 
U(1)' D symmetry is spontaneously broken in the weak coupling region j9 > 1. This means 
that we cannot study the continuum limit of the large N gauge theory by the EK model 
when D > 2, since we have to send (3 to infinity when we take the continuum limit. However, 
it is found that when D is even, one can twist the boundary condition of the EK model so 
that the U(1) D symmetry is not broken in the weak coupling region, while keeping the proof 



of the equivalence valid. This is the twisted Eguchi-Kawai (TEK) model [[T4j] defined by the 
action 

D 

S TEK = -N(3 Z, lu tr{U^U u UlUt}, (2.6) 

where is called as the 'twist', and defined as 

Z M „ = Z*^ = exp (2m/ L) for /i < u, (2.7) 

where N = L D I 2 . The observable that corresponds to the Wilson loop W(C) is now given 
by 

w{C) = ^(^Z p ^tv{U a U p U,---U x }, (2.8) 

where P^ u is the number of plaquettes in the (/ii/) direction on the surface spanned by the 
Wilson loop C. The TEK model has been recently used for searching for nontrivial fixed 
points in the six- dimensional large N gauge theory . 

In two dimensions, since the U(l) 2 symmetry is not spontaneously broken even in the 
weak coupling region, the EK model is equivalent to the lattice gauge theory in the planer 
limit in the above sense. The two-dimensional lattice gauge theory is solvable ]16|] in the 
planar limit and indeed the exact results obtained there has been reproduced by the EK 
model in the planar limit []nj. On the other hand, the TEK model can be considered in two 
dimensions as well, and its planar limit should be the same as that of the EK model. Note, 
however, that the two models can show different finite N effects, which is relevant to us, 
since we are aiming at discovering a large N limit other than the planar limit. Indeed we 
see that the approach to the planar limit from finite ./V is quite different for the two models, 
and the double scaling limit seems to exist only for the EK model. 



3 



Let us explain the algorithm we use for the Monte Carlo simulation of the EK and TEK 
models. Since the actions of the two models ( |2.1| ) and fl2.6| ) are not linear in terms of link 
variables unlike that of the ordinary lattice gauge theory ( |2.4| ), we cannot apply the heat 
bath method ]l7j as it stands. We therefore employ the technique proposed by Ref. [ |18|| . 
and introduce the auxiliary field Q^ v [p, < v), which are general complex N x N matrices, 
with the following action. 



S = Np^tiQ^Qr 

\1<V 

-N0 ]T trQt„ {t, v U,U v + t VIM U v U, 



fl<U 

-N(3 £ trQ M , (t; u UlUl + KJJttf) , (2.9) 



where t^ u is 1 for the EK model and is JZ^ V for the TEK model. Integrating out the 
auxiliary field Q^ U) one reproduces the original actions of the reduced models. The update 
of Q^v can be done easily by generating Gaussian variables. Note that the action ( |2.9| ) is 
linear in terms of U^, which means that we can apply the heat bath method. The update of 
Up is carried out by successively multiplying it by matrices each belonging to the N(N — 1)/2 
SU(2) subgroups of SU(iV) |L9||. We define 'one sweep' in this paper by the update of all 



the elements of Q^s and the update of all the U^s by multiplying elements of every SU(2) 
subgroup once for each. 

3 Planar limit and the finite N effects 

In this section, we see how the EK model and the TEK model at finite N approach the 
known results in the planar limit as one increases N with fixed (3 in the action ( [2.1|) and 
( |2.6|) , respectively. 

The observables calculated by Gross- Witten [Tj| for the two-dimensional SU(iV) lattice 
gauge theory in the planar limit is the Wilson loop W(C), where the contour C is a rectangle. 
Due to the discrete translational and rotational invariance and the parity invariance of the 
lattice gauge theory, the expectation value of Wilson loops is real and it depends only on its 
shape but not on how the loop is placed on the lattice nor on the orientation of the loop. 
We therefore denote the expectation value of the I x J Wilson loop by (W(I x J)). 

We examine what we get for the corresponding observables ( |2.5|) and (|2.8|) in the re- 
duced models with finite N. Just as we mentioned above for the lattice gauge theory, the 
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expectation value of the observable in the EK model that corresponds the Wilson loop is 
real and it depends only on its shape even for finite N due to the symmetries : U\ — > 
and U2 — > Ul, which corresponds to the rotational symmetry, and U\ — > U% and XJ% — > U\, 
which corresponds to the parity invariance. This is not the case, however, with the TEK 
model, in which the second one ('parity transformation') is not a symmetry. Accordingly, 
the expectation value of the observable in the finite N TEK model that corresponds to the 
Wilson loop is not real and it becomes complex conjugate when one flips the orientation 
of the loop, though the imaginary part should vanish in the planar limit, where it should 
reproduce the Wilson loop in the lattice gauge theory. When we define (w(C)) in reduced 
models which corresponds to (W(C)) in the lattice gauge theory, the real part is implicitly 
taken in the case of the TEK model. When we measure (w(C)) in either model, we average 
over the transformation of the Wilson loop by rotation and parity in order to increase the 
statistics. 

The internal energy of the system can be given by the expectation value of the lxl 
Wilson loop 

E(8) = {W(lxl)), (3.1) 



and the analytical result in the planar limit is obtained as |L6 



*<0 = K-A ^ 



P for 8<\ 

— for 8 > - 

4/3 iui r — 2" 

As one sees from this result, the second derivative of the internal energy is discontinuous at 
8 = 1/2, which means that the system undergoes a third-order phase transition at (5 = 1/2. 

In Fig. [I], we plot the internal energy against 8 for the two models for iV = 32. We 
can see that the data agree nicely to the theoretical result in the planar limit ( |3.2| ), which 
is represented by the solid line. In Ref. ||14|| , the results for N = 20 are given, which shows 
that the data for the EK model are slightly larger than the planar result, while the data 
for the TEK model are in complete agreement with it. This tendency has been observed in 
our simulation as well with smaller N. Thus as far as the internal energy is concerned, the 
convergence to the planar limit is better for the TEK model than for the EK model. 

Let us turn to the expectation value of rectangular Wilson loops of various size. The 
planar result is given by 

(W(I x J)) = exp (-k(8)IJ) , (3.3) 
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N=32 



1.0 



W 0.5 




Gross-Witten 
EK model 
□ TEK model 



0.5 

P 



Figure 1: The internal energy is plotted against (3 for the EK model (circles) and for the 
TEK model (squares) for N = 32. The solid line represents the theoretical result in the 
planar limit (N — > oo with fixed (3) obtained by Gross-Witten. 
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rectangular Wilson loop of various size 
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Figure 2: The expectation value of w(I x J) in the EK model (N — 16, (5 — 1.0) which 
corresponds to / x J Wilson loops is plotted against the physical area S = a 2 IJ. Each 
symbol represents the data for 1 = 1 (circles), 1 = 2 (squares), 1 = 3 (diamonds), 1 = 4 
(triangles) and 1 = 5 (tilted triangles) with J > I. The straight solid line represents the 
planar result by Gross- Witten. 
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where 

f "log/3 for f}<\ 

»W = {-log(l-4) fo r/ 3>l. < 3 ' 4 ' 

For I = J = 1, the above result reduces to the one for the internal energy ( |3.2| ). This 

result shows that the rectangular Wilson loops obey the area law exactly for all (3. From 

this result, one can figure out how to fine-tune the coupling constant (3 as a function of the 

lattice spacing a when one takes the continuum limit a — > 0. Since the physical area is given 

by S = a 2 IJ, we have to fine-tune (3 so that n([3)/a 2 is kept fixed. We therefore take 



/c(/5) (3.5) 



in the following. Since k((3) is given by eq. (|3.4|) , we have to send (3 to infinity as we take 
the a — > limit. 

In Fig. we plot the expectation value of I x J Wilson loops in the EK model against 
the physical area S = a 2 IJ for N = 1Q and (3 = 1.0, where I and J run from 1 to 16. 
Although the planar result depends only on the area of the Wilson loop, but not on its 
shape, we see that the results for finite N depend on the shape as well. For example, we see 
that the result for the 1x9 Wilson loop is farther from the planar limit than that for the 
3x3 Wilson loop. In general, one sees that the finite N effect for Wilson loops with the 
same area becomes larger as it becomes longer in one direction. This means that we have to 
specify not only the area but also the shape of the Wilson loop when we discuss the scaling 
behavior in the double scaling limit, as well as when we discuss the finite N effects. 

In Fig. |^ we plot the expectation value of square Wilson loops in the EK model against 
the physical area S = (al) 2 for (3 = 4.0 with N = 16, 32, 64 and 128. One can see that 
the data points approach the planar result from above monotonously. In Fig. f| we plot the 
deviation of the data for fixed Wilson loops from the planar result as a function of N for 
(3 = 1.0. One can see a clear 1/N 2 behavior, which means that the deviation is dominated 
by the subleading term in the large N limit. For (3 = 4.0, the deviation decreases slower 
than 1/N 2 up to iV = 128, which means that the contribution from the higher order terms 
is still comparable with the subleading term. The double scaling limit we consider in the 
next section corresponds to taking the large N limit together with the (3 — > oo limit so that 
all the terms in the 1/N 2 expansion contribute. 

The finite N effects for the TEK model appear in quite a different manner. Fig. |5] shows 
the results for the TEK model for the same set of values of N and with the same (3 as in Fig. 
|3l. The data point for the largest S for N = 128 is omitted, since the result turned out to be 



square Wilson loop in EK model 




Figure 3: The expectation value of w(I x I) in the EK model for /3 = 4.0 is plotted against 
the physical area S = (al) 2 . Each symbol represents the data for N = 16 (circles), N = 32 
(squares), N = 64 (diamonds) and iV = 128 (triangles). The straight solid line represents 
the planar result by Gross- Witten. 
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convergence in large N 

(3=1.0,EK model 

1 p _ ^ ^ ^ ^ — ^ — ^ 




10 100 
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Figure 4: The deviation of the expectation value of w(I x 7) in the EK model from the 
planar result is plotted against N for 7 = 1 (circles), 7 = 2 (squares) and 7 = 3 (diamonds) 
at P — 1.0. The solid straight line represents a fit to const. /N 2 , which is the dominant 
subleading term in the large N limit. 
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negative within the error bar. We have checked the thermalization by using hot starts and 
cold starts. The number of sweeps needed to achieve the thermalization is less than 1000 
sweeps for iV = 128. The errorbars are estimated using the standard jackknife method. The 
autocorrelation time is approximately 200 sweeps for N = 128. One sees from the results for 
fixed area with S < 2.5 that the data do not approach the planar result monotonously, but 
rather show some overshooting behavior. For larger S, the tendency of convergence is not 
yet seen up to N = 128. Thus, although it might be that the convergence to the planar limit 
is somewhat accelerated in the TEK model as is seen from the results for smaller Wilson 
loops, we find that the finite N effects in the TEK model are much more complicated than 
in the EK model. 

Another interesting quantity known in the planar limit is the expectation value of Wilson 
loops with self-intersections | 20| . In Fig. |6], we illustrate a typical Wilson loop we consider. 



It is given by the trace of the product of link variables along the loop C = abcdaefghia. We 
split the original loop C into two loops at the self-intersection point and define C\ = abcda 
and C 2 = aefghia. 

The expectation value of such a Wilson loop has been calculated by solving the 



Makeenko-Migdal equation [21] in the continuum and the result is given by 



(W{C)) = (1 - 2KS{C X )) exp [-K{S{Cx) + S{C 2 )}} . (3.6) 

S(C) is the physical area of the region surrounded by the loop C. When S{C\) = 0, the 
above expression reduces to the one for the Wilson loop without self-interactions, which 
obeys the area law. K is the physical string tension, which is taken to be 1 in this paper, 
since we have fixed the scale by eq. ( |3.5|) . When S(Ci) > 0, the factor (1 — 2KS{C\)) gives 
a nontrivial correction. We measure the observable 

w(C) = Itr (u^Utu\ L Ul Ll U^Ul L Ul L2 ) , L, < L 2 , (3.7) 

in the EK model and the corresponding one in the TEK model, where we multiply the above 
expression by Z^ 2 Ll+L ^ according to the definition ( |2.8| ). Comparison with the continuum 
result in the planar limit ( |3.6|) can be made by putting S(C\) = a 2 LLi and S(C2) = a 2 LL 2 . 
In order to extract the nontrivial factor (1 — 2KS(C\)), where K — 1, we measure the expec- 
tation value (u>(C)) with fixed S{C\) + S{C 2 ) and plot the result against S{C\). Specifically, 
we take L = 1 and consider Li and L 2 with L\ + L 2 = 8 and L\ = 0, 1, 2, 3, 4. The observable 
for L\ = is nothing but the (L x L 2 ) Wilson loop. Fig. |7| and Fig. [8] shows the results 
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square Wilson loop in TEK model 




Figure 5: The expectation value of w(I x J) in the TEK model for (3 = 4.0 is plotted against 
the physical area S = (al) 2 . Each symbol represents the data for N = 16 (circles), N = 32 
(squares), N = 64 (diamonds) and iV = 128 (triangles). The straight solid line represents 
the planar result by Gross- Witten. 
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Figure 6: A typical Wilson loop with a self- interaction 

for N = 32 and N = 64 respectively. The planar results by Kazakov-Kostov are shown by 
the straight heavy solid lines. One sees that the data are on the lines parallel to the planar 
results, and the lines are slightly above the planar results for the EK model, and slightly 
below the planar results for the TEK model []. The discrepancy is smaller for N = 64 than 
for N = 32. 

This shows that the finite N effect comes mainly from the area law part, and the correc- 
tion due to the factor (1 — 2KS(Ci)) in (|3.6| ) does not show any significant finite N effect. 
We also note that the total physical area S(Ci) + S(C2) is 0.516 for (3 = 4.0 and the tendency 
of convergence to the planar limit is consistent with that observed in Fig. |3] and Fig. |5|. 

4 Existence of double scaling limit in EK model 

In this section, we consider the EK model not as a reduced model which reproduce the lattice 

gauge theory in the planar limit, but as a matrix model which defines a string theory. We 

therefore examine whether we can take a large N limit other than the planar limit discussed 

in the previous section. 

We consider the EK model as a toy model of the IIB matrix model. Just as in Ref . ||, |l(| , 

we make the T-duality transformation, when we interpret the EK model as a string theory. 

Upon the T-duality transformation, the compactification radius of the torus to which the 

space time is compactified is inversed and the winding mode and the momentum mode are 

exchanged. Since the EK model can be considered to be defined on a unit cell of size a with 

the periodic boundary condition, it can be considered, after the T-duality transformation, 

Strictly speaking, L, L\ and Li should be large enough when we compare the data with the continuum 
result. The data show, however, that the finite lattice spacing effect is absent accidentally for the observable 
considered, as in the case with the simple rectangular Wilson loops, which obey the area law exactly even 
for the (1 x 1) loop. 
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Wilson loop with self-intersection 

N=32 




Figure 7: The expectation value of w(C) in the EK (open symbols) and the TEK model 
(filled symbols) for N = 32 with (3 = 4.0 (diamonds), 6.0 (squares), 8.0 (triangles) and 10.0 
(circles) are plotted against the physical area S(Ci) = a 2 LLi. The straight heavy solid lines 
represent the continuum results in the planar limit by Kazakov-Kostov. 
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Wilson loop with self-intersection 

N=64 




Figure 8: The expectation value of w(C) in the EK (open symbols) and the TEK model 
(filled symbols) for N = 64 with (3 = 4.0 (diamonds), 6.0 (squares), 8.0 (triangles) and 10.0 
(circles) are plotted against the physical area S(Ci) = a 2 LLi. The straight heavy solid lines 
represent the continuum results in the planar limit by Kazakov-Kostov. 
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as a string theory in the two-dimensional space time compactified on a torus of size 1/a. 
The observables we consider in this section is the Wilson loops 



tr {U a UpUT--Ux}, (4.1) 

where the suffix a, /?,■••, A runs over ±1,±2 and E/_ a (a > 0) is defined by C/_ Q = C/t. The 
winding number of the Wilson loop is given by = d M + /3 M + • • • + A M , where d M denotes 
a unit vector in the a direction. Note that the observables ( |2.5|) considered in the previous 
section can be viewed as the Wilson loops in the EK model with = 0, while in this section 
we consider those with n M ^ as well. After the T-duality transformation, the winding of 
the Wilson loops represents the momentum distribution of the string in the two-dimensional 
space time. The total physical momentum carried by the string which corresponds to the 
Wilson loop (|4.1| ) is given by P M = n^a. Under the U(l) 2 transformation ( |2.2|) , the Wilson 



loop (4.1) is multiplied by exp(m^6^). Since the symmetry is not spontaneously broken 
in the present model, the one-point function of Wilson loops in the winding sector should 
vanish. Similarly, the n-point function should vanish in general, unless the sum of the 
winding number of the Wilson loops is zero. After the T-duality transformation, the above 
selection rule gives nothing but the momentum conservation in the two-dimensional space 
time. This corresponds to the fact that the U(l) 2 symmetry is nothing but the translational 
invariance when we identify the phase of the eigenvalues of as the two-dimensional space- 
time coordinates. 

In what follows, we consider 

{ $(L) = tr(U^Ut a U^) 

OW(L) = tr(t£), (4.2) 

which give typical Wilson loops in the non-winding sector and the winding sector respectively. 
Due to the trace, the independent operators in the non-winding sector are actually given by 
0^2 {L) and 0%l\L), which are complex conjugate to each other. When we take the a — > 
limit, we have to send L to the infinity by fixing aL. Note that in ( |4. 1| ) and fl4.2p , we do 
not put 1/JV in front of the trace as we did in (|2.5| ). This is because with this convention, 
the power counting of N for correlation functions matches with the one conventional in the 
old-fashioned matrix model |l[|. 

We do not consider the TEK model in what follows, since the interpretation as a string 
theory is not as obvious as in the EK model. Indeed, as is seen in Fig. |5], the TEK model 
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does not show any systematic finite N effects, while the EK model does. It is therefore 
natural to consider that the TEK model does not allow any sensible double scaling limit. 

Let us first examine the one-point function of Wilson loops in the non-winding sector, 
namely Gi(L) = (0^(L)). The observable (w(L x L)) considered in the previous section 
is nothing but N~ 1 G ( ( l \L). We have seen in Fig. ^ how the data for iV _1 G^(L) approach 
the planar limit when we increase N with fixed (3. In Fig. |9|, we plot N~ 1 G ( ( l \L) against 
S = {aL) 2 for fixed iV with various values of (3. One can see that the effect of increasing (3 
is quite similar to the finite N effect. Thus it is natural to expect that we may be able to 
balance these two effects by taking the limit (3 — > oo together with the N — > oo limit. Indeed 
we find that N~ l G^\L) scales in the large N limit by fixing (3/N . In Fig. [K| we show the 
data for (3/N = 1.0/32 and (3/N = 1.5/32. N is taken to be 32, 48, 64, 96 and 128. We see a 
clear scaling behavior. Measurements have been done every 10 sweeps and we average over 
2000 configurations for iV = 32 ~ 96, and over 1000 configurations for N = 128. All the 
data in this section have been obtained in this way, unless mentioned. 

The finite N effect is seen as a deviation from the scaling behavior for large S. The 
deviation shows up at L ~ N/8. Since a oc l/>/j3~ oc 1/yN for large N, the scaling region 
enlarges as \fN . The scaling function for fixed (3/N approaches the planar limit as (3/N 
becomes smaller. 

From this observation, it is natural to expect that the double scaling limit of the EK 
model can be taken by fixing (3/N, which corresponds to the string coupling constant g str . 
The planar limit corresponds to g str = 0. 

Let us next consider Wilson loops in the winding sector. Since the one-point function 
vanishes as is explained above, we consider the two-point function 

dF\L) = (Ol?\L)0^(L)). (4.3) 

Note that the right hand side does not depend on /i, due to the rotational invariance. In the 
measurement, we average over /i = 1,2 to increase the statistics. This quantity is naively 
0(N 2 ), but due to the factorization, its leading term is given by (0^(L))(Oi^(L)), which 
actually vanishes due to the U(l) 2 symmetry. The order of the quantity is therefore 0(1). 
We can regard the 0(1) quantity as the connected two-point function of the Wilson loops 
in the winding sector. 

Let us examine if the two-point function G2 (L) allows the double scaling limit. In Fig. 
II], we plot G { 2 ] {L) against I = aL for (3/N = 1.0/32 with N = 32,48,64,96 and 128. No 
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one-point function in the non- winding sector 



N=32, (3=1.0,2.0,3.0,4.0 




Figure 9: The one-point function of the Wilson loops in the non-winding sector N~ 1 G 1 (L) 
in the EK model is plotted against S = (aL) 2 for N = 32 with f3 = 1.0 (circles), 2.0 (squares), 
3.0 (diamonds) and 4.0 (triangles). The straight solid line represents the planar result by 
Gross- Witten. 
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one-point function in the non- winding sector 
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Figure 10: The one-point function of the Wilson loops in the non- winding sector N~ 1 G{ (L) 
in the EK model is plotted against S = {aL) 2 for the sets of f3 and N with f3/N fixed to 
1.0/32 and 1.5/32. Each symbol represents the data for iV = 32 (circles), N = 48 (squares), 
iV = 64 (diamonds), N = 96 (triangles) and iV = 128 (inverse triangles). The straight solid 
line represents the planar result by Gross- Witten. The dashed lines, which connect the data 
for iV = 128, are drawn to guide the eye. 
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scaling behavior is observed. However, as is the case with the double scaling limit in the 
old-fashioned matrix models JTJ, we might have to admit the wave function renormalization. 
Since the Fig. [Tl| is a log-log plot, the wave function renormalization amounts to shifting 



the curves relatively along the vertical axis. Fig. [12] shows the result for the renormalized 
two-point function 

G^\L)=p- 0S5 G^\L). (4.4) 

One can see a clear scaling behavior in the intermediate region of I. The scaling function 
is a power of / for large I, and the power is 0.80(2). The finite N effect can be seen as a 
saturation of the power-law behavior in the large I region, which shows up at L ~ N. This 
is to be compared with the case with iV _1 G^(L) in Fig. [10], where the deviation from the 
scaling behavior in the large / region shows up at L ~ N/8. This discrepancy is not strange 
since the meaning of L depends on which of the two operators O^(L) and 0^\V) one 
considers. 

Unlike in Fig. pi], there is a deviation from the scaling behavior also in the small / region. 
The tendency of convergence to a scaling function in this region is not clearly seen up to 
iV = 128. The fact that such a deviation does not appear for A^~ 1 G^(L) in Fig. [1^ might 
be rather considered as an accidental property of the Wilson loop in the non-winding sector. 

We can see the scaling behavior for different j3/N as well. In Fig. [T3], we plot the renor- 
malized two-point function G^ } (L) against / = aL for (3/N = 1.5/32 with N = 32, 48, 64, 96 
and 128. The behavior is qualitatively the same as that of Fig. [12|. Note that this result 
shows that the wave function renormalization does not depend on g str = {3/N, which is to 
be expected. The power of the scaling function in the large I region is 0.72(3). 

In order to understand the power- law behavior, we calculate G^"\l) for (3 = in the 
large N limit analytically through Schwinger- Dyson equation Q. We consider the quantity 
(tr(t a U^)tr(Uj 1 L )) , where t a (a — 1, ■ • • , N 2 ) denote the generators of the U(iV) group. By 
changing the variable of integration as — > (1 — iet a )U fl , we obtain the following identity. 

£ (tr(t a U^ M t a U^)tr(Ul L )) - L(tr(t a U^tT(t a Ul L )} = (4.5) 

M=l 

Taking a summation over a = 1, • • • , iV 2 , we obtain 

NMUZMU?)) + E (tr(^- M )tr(^ M )tr(^)) _ NL = , (4.6) 

M=l 



2 We thank H. Kawai for suggesting the use of Schwinger-Dyson equation for this purpose. 
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two-point function in the winding sector 

without wave function renormalization 
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Figure 11: The two-point function in the winding sector without the wave function renor- 
malization G^"\l) is plotted against I = aL for sets of f3 and N with (3/N = 1.0/32. Each 
symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 96 (triangles) and N = 128 (inverse triangles). 
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two-point function in the winding sector 

with wave function renormalization 
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Figure 12: The two-point function in the winding sector with the wave function renormal- 
ization G^"\l) is plotted against I = aL for sets of (3 and N with f3/N = 1.0/32. Each 
symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 9Q (triangles) and N = 128 (inverse triangles). 
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two-point function in the winding sector 

with wave function renormalization 
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Figure 13: The two-point function in the winding sector with the wave function renormal- 
ization Gq"\l) is plotted against I = aL for sets of (3 and N with (3/N = 1.5/32. Each 
symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 96 (triangles) and iV = 128 (inverse triangles). The result for (3/N = 1.0/32 with 
iV = 128 are replotted by the heavy line. The dotted line is a line with slope 1, which 
represents the result for (3 = and N = oo. 
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where we have used the identity Y.a^ a )ij{^ a )ki = SuSjk- Since the three-point function 
(tT(U^~ M )tr(U^ I )tr(Uj l L )) is at most 0(1/N) as we explain later, we obtain 

Gj 0) (L) = (tr(E^)tr(^)> = L. (4.7) 

We have checked this result explicitly by Monte Carlo simulation. It is natural to expect 
that in the planar limit the two-point function in the large I region can be described by the 
strong coupling limit obtained above. The fact that the power of I in the scaling function 
approaches one for decreasing g str = {3/N can be naturally understood in this way. 

Note that the wave function renormalization was actually necessary for G { r\L) as well, 
since we had to multiply it by iV _1 when we see the scaling behavior. In this sense, we 
should say that the renormalized one-point function is given by Gi (L) = /3 -1 Gj n (L). 
Thus, although the result for the two-point function of the Wilson loops in the winding 
sector supports to some extent our conjecture concerning the existence of the double scaling 
limit, we encounter a new problem: why the wave function renormalization for the one-point 
function G^lL) is not equal to that for the two-point function G%{Ij). 

In order to clarify this problem, we consider the two-point function of Wilson loops in the 
non-winding sector defined by {0$(L)0$(L)) and (0$ (L)Oz± (L)) . These observables are 
0(N 2 ) quantity, but the leading term is given by the disconnected part (o{f '(L)) (o[f (L)) 
and (O12 (L)) (O21 (L)} , respectively, due to the factorization. The connected two-point 
function can be defined by 

(o£>(L)0®(L)) c = (0^(L)0S>(L))-{0^(L)){0S(L)) (4.8) 
(0£>(L)0%\L)) c = (0%\l)0£\l))-(0&\l))(0£ ) (L)), (4.9) 

which are 0(1) quantities. 

Since (O^(L)} = (0%l\L)) G TZ, as is mentioned in the previous section, we can avoid 
the subtraction by considering the difference of the two-point functions. 

G^-(L) = ±{(0%\l)0%\l)) c -(0%\l)0%\l)) c } (4.10) 

= ~ {{0^(L)0^(L)) - (0®{L)0®(L))} (4.11) 
= ((ImO^(L)) 2 ). (4.12) 

We have used the fact that (0$(L)0$(L)) = (0{f (L)0$> (L)) G K in the third equality. 
In Fig. M we plot the G^ n)_ (L) against S = (aL) 2 for /3/N = 1.0/32 with N = 32, 48, 64, 96 
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two-point function in the non- winding sector 

without wave function renormalization 
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Figure 14: The two-point function in the non-winding sector without the wave function 
renormalization G^ - (L) is plotted against S = (aL) 2 for sets of (3 and iV with (3/N = 
1.0/32. Each symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 
(diamonds), N = 96 (triangles) and iV = 128 (inverse triangles). The solid line, which 
connects the data for iV = 128, is drawn to guide the eye. 
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two-point function in the non- winding sector 

with wave function renormalization 
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Figure 15: The two-point function in the non-winding sector with the wave function renor- 
malization G^~{L) is plotted against S = (aL) 2 for sets of (3 and iV with (3/N = 1.0/32. 
Each symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 96 (triangles) and N = 128 (inverse triangles). The solid line, which connects the data 
for iV = 128, is drawn to guide the eye. 
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and 128. There seems to be no scaling behavior. Let us assume the same wave function renor- 
malization as the one for G^\l) in (O). In Fig. pL we plot G^' (L) = (3~ om G^ {L) 



for p/N = 1.0/32 with JV = 32,48,64,96 and 128. We find a clear scaling behavior. The 
discrepancy in the large S region starts from L ~ N/8, which is the same as with G^iL). 
This suggests that the operator O^(L) gives a finite A" effect to the correlation function 
which includes it for L ~ N/8 in general. It is natural to expect that this is the case also 
with the operator 0^"\L) with L ~ N, which we exploit in discussing the scaling behavior 
of the three-point function of Wilson loops in the winding sector later. The scaling seems to 
extend to the small / region except for the data that corresponds to L = 1. This is consistent 
with our previous observation that the Wilson loops in the non-winding sector do not suffer 
from finite A" effects in the small I region, unlike those in the winding sector. 
Similarly, we can define 

Gf )+ {L) = \{(0^\L)0^\L)) c+ (0^\L)0^(L)} c } (4.13) 
= ((ReO^L)) 2 )-^^)) 2 . (4.14) 

The result for G^/^iL) = /J? -0,65 ^™^ (L) also shows the scaling behavior, which is qualita- 
tively the same as in Fig. [T5[ albeit larger error bars due to the subtraction. 

As in the case of G^"\l), we can calculate G% n (L) for (3 = in the large A" limit 
by using the Schwinger- Dyson equation and obtain G2 (L) = 1/2. This suggests that in 
the planar limit, the asymptotic behavior of the scaling function of G% ™ (L) for large I is 
constant, which is consistent with our results. 

From these observations, we conclude that only the one-point function is exceptional 
concerning the wave function renormalization. Note that the wave function renormalization 
of GS n) (L) is, in some sense, kinematically constrained by the fact that G^iO) = N. The 
non-vanishing one-point function can be considered as the background of the corresponding 
string field and the exceptional behavior is not very unnatural. 

As a further check of the existence of the double scaling limit, we consider the three-point 
function of Wilson loops in the winding sector defined by 

G$°\L) = (0^(L)0^(L)0^(2L)}. (4.15) 

Note that the right hand side does not depend on /i, and we average over ji = ±1,±2 in 
the measurement to order to increase the statistics. The three-point function is naively 
C(A^ 3 ) quantity, but the leading term (O^ (L)) (Oj°\L)) (0 {w) J2L)) and the subleading 
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three-point function in the winding sector 

without wave function renormalization 
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Figure 16: The three-point function in the winding sector without the wave function 
renormalization G^°\l) is plotted against / = aL for sets of (3 and N with (3/N = 1.0/32. 
Each symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 9Q (triangles) and N = 128 (inverse triangles). 
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three-point function in the winding sector 



with wave function renormalization 
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Figure 17: The three-point function in the winding sector with the wave function renor- 
malization G^\l) is plotted against I = aL for sets of f3 and N with (3/N = 1.0/32. Each 
symbol represents the data for N = 32 (circles), N = 48 (squares), N = 64 (diamonds), 
N = 96 (triangles) and iV = 128 (inverse triangles). 
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term {0^\L)0^\L)) (CrJ(2L)) etc., which are of order 0(N), vanish due to the U(l) 2 
symmetry, and it is actually 0(1/ N) quantity, which can be regarded as the connected 
three-point function. 



In Fig. |T|and Fig. we plot Gf ] (L) and G^(L) = (/?/4.0)- a65x3 / 2 G^ ] (L), respec- 



tively, against I = aL for (3/N = 1.0/32 with N = 32,48,64,96 and 128. Measurements 
have been done every 10 sweeps. Each point is an average over 4000 configurations for 
iV = 32 ~ 96, and over 2000 configurations for N = 128. Note that since the data can be 
both positive and negative, we cannot make a log plot, though we should do so when we dis- 
cuss scaling behaviors. In order to avoid a possible unfair comparison of the figures, we have 
defined the G3 (X) so that the data for iV = 128 look the same as those for G { ™\l). Al- 
though the data are too noisy to confirm the scaling behavior exclusively with the particular 
wave function renormalization assumed here, the renormalized data are at least consistent 
with the scaling behavior in some region of /. Furthermore, the deviation from the scaling 
in the large I region shows up at 2L ~ N, which is in agreement with the finite N effects 
seen for the two-point function of the same operator, namely G 2 W (L). This is not the case 
with the data without the wave function renormalization. We therefore conclude that the 
data for the three-point function also support the existence of the double scaling limit with 
the universal wave function renormalization. 

5 Conclusion and discussion 

In this paper, we studied the EK model as a toy model of the IIB matrix model. In the 
planar limit, the one-point function of the Wilson loops in the no n- winding sector gives the 
expectation value of Wilson loops in the large N lattice gauge theory in two dimensions. 
Using the exact results in the 2D large N lattice gauge theory, we fixed how to send a to 
as we send (3 to infinity as 



In the planar limit, N is sent to infinity before taking the limit a — > and (3 —>■ 00. We 
found that there is another way of taking the large N limit, namely with f}/N fixed. We 
found that the correlation functions of Wilson loops in the non-winding sector as well as in 
the winding sector have nontrivial limits with the wave function renormalization / ^ -0 - 65 / 2 for 
each Wilson loop. The only exception we saw is the one-point function, for which the wave 




(5.1) 
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function renormalization is (3~ l . We consider that this is not so unnatural, since the one- 
point function is special in the sense it can be considered as the non-vanishing background 
of the corresponding string field. 

The deviation from the scaling behavior due to the finite N is seen in the large / region. 
The number of the links for which the Wilson loop operator gives finite N effects depends 
on the operator considered. It is around N for the Wilson loops in the winding sector and 
around N/2 for the Wilson loops in the non- winding sector. In either case, the physical 
scaling region enlarges as y/N in the large / region. 

On the other hand, the situation in the small I region is not so clear. The deviation 
from the scaling behavior is seen for the two-point functions of Wilson loops in the winding 
sector. Convergence to a scaling function seems to be slow. For the correlation functions of 
the Wilson loops in the non-winding sector, on the other hand, no significant deviation from 
the scaling behavior is seen in the small / region. We consider this as an accidental property 
of Wilson loops in the non-winding sector. 

The data for the n-point functions with n > 3 are unfortunately too noisy to confirm 
the scaling behavior for fixed f3/N. We found, however, that the data for the three-point 
function of Wilson loops in the winding sector with the same wave function renormalization 
are consistent with the scaling. Furthermore, the deviation from the scaling in the large / 
region due to the finite N effects appears in a manner which is consistent with the finite N 
effects seen for the two-point function of the same operator. 

We calculated the two-point functions for (3 = in the large N limit by using the 
Schwinger-Dyson equation. The result seems to describe the large I behavior of the corre- 
sponding quantity in the planar limit. We can calculate the connected n-point functions in 
general, which are order 0(1/N^ n ~ 2 ^) quantities. We find that the coefficient of l/iV (n ~ 2 ) is 
zero for j3 = in the large N limit. This suggests that the n-point functions with n > 3 
in the planar limit goes to zero in the large / region. One might fear that the string theory 
constructed through the double scaling limit of the EK model is a free theory. Our data 
suggest, however, that the three-point function is not constantly zero at least when g str ^ 0. 

A technically important comment is that while the old-fashioned matrix model, whose 
action is unbounded from below, becomes well-defined only in the large N limit, the EK 
model as well as the IIB matrix model is completely well-defined even for a finite N. This 
makes these models much easier to study numerically than the old-fashioned matrix model. 

To summarize, our numerical results suggest strongly that the double scaling limit of 
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the EK model can be taken by sending iV to the infinity with the following combinations of 
parameters fixed. 

P/N ~ g str (5.2) 
Na 2 ~ (a'g str )-\ (5.3) 

where g s t r is the string coupling constant and a' is the string tension. 

It is intriguing to ask what is the string theory constructed through the double scaling 
limit of the EK model. We interpreted the model as a string theory through the T-duality 
transformation as is done in the IIB matrix model. The target space is therefore two- 
dimensional space time compactified to a torus of size ~ 1/a, which goes to infinity as 
a — > 0. The theory must be different from the c = 1 matrix model interpreted as 2D 
string theory with linear dilaton background, since we have at least the discrete rotational 
invariance in the EK model. Note also that although the target space is two-dimensional, 
the dynamical degrees of freedom of the link variables are actually absorbed by the gauge 
degrees of freedom. It might be that this theory cannot be described as 2D gravity coupled 
to some conformal matter. 

An important difference between the EK model and the IIB matrix model is that the 
dynamical variables in the former are the unitary matrices, while those in the latter are the 
hermitian matrices. Considering the interpretation of the models as string theories using 
the T-duality transformation, it seems to be more natural to formulate the matrix models 
in terms of unitary matrices. Constructing such a formulation for type IIB superstring has 
been considered in Ref. |]22|| . The problem corresponding to the fermion doublers in the 
lattice chiral gauge theory has been overcome by the use of the overlap formalism, but the 
remaining problem is whether one can recover the supersymmetry in the large N limit with 
the proposed model or with fine-tuning of the coefficients of some possible counterterms. 
On the other hand, one could consider a hermitian matrix version of the two-dimensional 
EK model, with the action S = — ^trLA^, A u } 2 , but the theory is not well defined PBfl . We 
can make the theory well defined, for example, by adding a term like m 2 (tYA 2 l — -^(trA^) 2 ), 
which preserves the U(l) 2 symmetry : A^ — > A^ + a^. However, the U(l) 2 symmetry is 
then spontaneously broken and the model must be in a universality class other than the one 
corresponding to the two-dimensional EK model. 

Despite the above subtlety, the relations ( |5.2|) and ( |5.3| ) which specify how one should 
take the double scaling limit in the two-dimensional EK model might be naively compared 
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with the ones conjectured for the IIB matrix model [ID]; namely, 

g 2 N ~ a' 2 (5.4) 
Ne 2 ~ (5.5) 

e is a dimensionless constant which corresponds naively to a. We therefore have a ~ ea' -1 / 2 . 
This means that the relation (|5.5|) coincides with the corresponding one in the EK model, 
namely (|5.3|) . The coupling constant g can be related to (3 naively in the following way. We 
relate the hermitian matrices in the IIB matrix model to the link variables in the 
EK model through = exp(mA At ), and expand the action ( |2.1| ) in terms of a. We obtain 
S = — Nj3aHi[A^, A u ] 2 ■ ■ ■. In Ref. ||1 Q|| , the factor in front of the trace is denoted as l/4g 2 , 
which means that 1/g 2 ~ Nj3a 4 . Therefore, (|5.4j ) corresponds naively to fixing j3/N 2 in our 
notation. 

We also note that the density of eigenvalues is constant in the double scaling limit. As we 
mentioned above, the extent of the space time is given R = — . Since we have N eigenvalues 
distributed in the two-dimensional space time with this extent, the average denstity is given 
by 

P = ^ = (2^' (5 ' 6) 
which is constant in the double scaling limit. This fact is natural from the string theoretical 
point of view, since it means that there are only finite dynamical degrees of freedom on the 
average in a finite region of the space time. 

The fact that a sensible double scaling limit can be taken for the large N reduced model of 
two-dimensional lattice gauge theory is itself encouraging for the research of nonperturbative 
formulation of superstring theory through the IIB matrix model. We hope to report on 
numerical studies of the IIB matrix model in future publications. 
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